A model for microinstability destabilization and 
enhanced transport in the presence of shielded 3-D 
magnetic perturbations 

T M Bird 1 , C C Hegna 2 

1 Max-Planck- Institut fur Plasmaphysik, D- 17491 Greifswald, Germany 

2 Departments of Engineering Physics and Physics, 
University of Wisconsin, Madison, WI 53706 

E-mail: tbird@ipp.mpg.de 

Abstract. A mechanism is presented that suggests shielded 3-D magnetic 
perturbations can destabilize microinstabilities and enhance the associated anomalous 
transport. Using local 3-D equilibrium theory, shaped tokamak equilibria with small 
3-D deformations are constructed. In the vicinity of rational magnetic surfaces, the 
infinitc-n ideal MHD ballooning stability boundary is strongly perturbed by the 3- 
D modulations of the local magnetic shear associated with the presence of near- 
resonant Pfirsch-Schluter currents. These currents are driven by 3-D components of 
the magnetic field spectrum even when there is no resonant radial component. The 
infinitc-n ideal ballooning stability boundary is often used as a proxy for the onset of 
virulent kinetic ballooning modes (KBM) and associated stiff transport. These results 
suggest that the achievable pressure gradient may be lowered in the vicinity of low 
order rational surfaces when 3-D magnetic perturbations are applied. This mechanism 
may provide an explanation for the observed reduction in the peak pressure gradient 
at the top of the edge pedestal during experiments where edge localized modes have 
been completely suppressed by applied 3-D magnetic fields. 
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1. Introduction 

The use of applied three-dimensional magnetic fields can have beneficial effects on the 
performance of otherwise axisymmetric toroidal confinement devices [H [2] . In particular, 
the application of external 3-D resonant magnetic perturbations (RMP) whose resonant 
surface lies in the edge region of H-mode tokamaks can suppress the appearance of 
edge localized modes (ELMs) under certain conditions [3j HJ Ej. ELMs repeatedly and 
violently eject hot plasma from the tokamak edge onto solid materials of the device. In 
order to safely operate at the desired parameters without prohibitive erosion of first wall 
materials, future experiments such as ITER will require a reliable method for controlling 
edge properties to prevent large type-I ELMs [6]. The original motivation for pursuing 
the RMP method of ELM suppression was based on the expectation that resonant 
magnetic perturbations of sufficient magnitude would produce overlapping magnetic 
islands that would greatly enhance plasma transport in the edge region and accordingly 
reduce the drive for ELM inducing instabilities associated with edge pressure gradients 
and currents [7J [5]. However, plasma flows in the edge region may be of sufficient 
magnitude to shield resonant magnetic perturbations from penetration (9j [10] in which 
case stochastic transport is no longer viable. Hence, a different mechanism may be 
required to explain the observations which indicate that the externally applied 3-D 
perturbations can enhance plasma transport in the pedestal and preclude the triggering 
of ELMs. 

Experiments with applied 3-D magnetic perturbations have now been performed on 
a number of tokamaks worldwide, yielding varying results. Experiments performed at 
the DIII-D tokamak with n = 3 perturbations have demonstrated suppression of large 
type-I ELMs at various collisionalities and with different plasma shaping. In many cases 
the large type-I ELMs are replaced by higher frequency, lower amplitude edge transport 
bursts [llj . This has also been achieved for example at JET [12] with n = 1 and n = 2 
perturbations, at MAST with n = 3 perburations [13], and at ASDEX-U with n = 2 
perburbations [H]. However, in some DIII-D experiments a complete elimination of 
bursty edge transport events has even been achieved [5]. In this work we provide a 
possible explanation for ELM-free operation that does not rely on stochastic magnetic 
fields to enhance the transport in the edge region of H-mode tokamaks. 

We propose a model for the enhanced transport in the edge region that is 
based on the three-dimensional distortions of the magnetic surfaces adversely affecting 
microinstabilities. While the effects of 3-D fields on stability discussed here is generic to 
a broad range of localized modes, we concentrate on the properties of kinetic ballooning 
modes (KBMs) in this work, which are primary candidates for explaining edge pedestal 
transport in H-mode plasmas. KBM growth rates become very large as the infinite-n 
ideal ballooning boundary is approached [16j [15j [TF] . As such, ideal MHD ballooning 
mode calculations are often used as a proxy for the abrupt onset of KBM instabilities and 
the associated stiff transport response [18] . In particular, calculations using information 
from ideal ballooning stability and peeling-ballooning (P-B) theory used to predict 
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ELM onset have successfully predicted self-consistent pedestal widths and heights [T9] . 
Gyrokinetic modeling of microinstabilities in the pedestal also supports the hypothesis 
that KBMs limit the pedestal pressure gradient and thus regulate the inward advance 
of the pedestal [20|l2T]. 

In this work, the effect of the 3-D distortions of the MHD equilibrium on ideal 
ballooning stability is calculated. For amplitudes of comparable measure to those present 
in experiments, the 3-D distortions are shown to substantially modify the local shear 
and reduce the critical pressure gradient for ideal ballooning instability. This is purely 
a 3-D MHD equilibrium effect, occuring due to Pfirsch-Schliiter currents near rational 
surfaces that are driven by 3-D inhomogeneity in the magnetic field spectrum. This 
is a resonant effect and would lower the achievable pressure gradient in the vicinity of 
low order rational surfaces where the 3-D deformation is strongest. Enhanced transport 
due to the onset of KBM instabilities at lower pressure gradient, due to the 3-D effect, 
could halt the inward progression of the edge pedestal and achieve ELM mitigation 
by preventing the pedestal height from reaching values where P-B modes are driven 
unstable. This explanation is consistent with recent pedestal modeling work where the 
prevention of ELMs was explained by the existence of a 'wall' at a low order rational 
surface which prevents the inward advance of the pedestal [22] . 

It is worth noting that the global MHD instabilities that are thought to trigger 
ELMs, Peeling-Ballooning modes, may be sensitive to this local shear modulation as 
well. The high-n components of these modes are dominated by ballooning structure, and 
we show in this work that (local) infinite-n ballooning modes are strongly affected by the 
local shear modulation (though we use the infinite-n ballooning calculations merely as a 
proxy for the onset of KBM instabilities and the associated transport). To understand 
the effect of resonant Pfirsch-Schliiter currents on the Peeling-Ballooning instabilities 
would require global MHD stability calculations using global 3-D MHD equilibria which 
include the resonant Pfirsch-Schliiter physics. This is beyond the scope of this work. In 
particular, this model attempts to provide an explanation for the DIII-D experiments 
where the increased pedestal transport completely precludes the triggering of any ELMs 
or other bursty events. Therefore, we focus on local KBM stability (with infinite-n MHD 
ballooning mode stability as a proxy) and utilize radially local equilibrium, which greatly 
simplifies the modeling effort and should be adequate for microinstability calculations. 

In the following section, details for how 3-D fields alter the MHD equilibrium are 
described using local 3-D equilibrium theory. In section III, the infinite-n ideal MHD 
ballooning stability boundary is evaluated for a set of equilibria with 3-D flux surface 
deformations of experimentally relevant magnitude. It is found that as the q value 
approaches a rational value, the instability boundary is strongly modified by the presence 
of the 3-D fields. Section IV explains this result via a detailed examination of the 
local magnetic shear. The stability boundary is expanded when near-resonant Pfirsch- 
Schliiter currents modulate the local magnetic shear in a manner strongly conducive to 
ballooning instability. Section V summarizes these results and discusses the implications 
for understanding experiments where RMP suppression of Edge Localized Modes has 
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been achieved. Appendix A elucidates the role of the 3-D fields on the local shear and 
provides an analytic estimate for the critical pressure gradient for ballooning instability 
for a low-/?, high aspect ratio, circular cross section tokamak when weak 3-D fields are 
present. 

2. The local 3-D equilibrium model 

A description of the equilibrium in the presence of 3-D fields can be provided by the 
prescriptions of local 3-D MHD equilibrium theory (23]. In this formulation, the shape 
of the 3-D magnetic surface is parametrized by straight-field line angles, the value of 
the safety factor q and two profile quantities, typically the presure gradient and the 
average magnetic shear q'. These parameters uniquely specify the properties of 3-D 
MHD equilibrium in the vicinity of a flux surface. The original motivation for this 
formulation was to study the effects of 3-D shaping on ballooning stability in stellarator 
configurations [23], and it has also proved useful for interpreting LHD data [25]. This 
formulation provides an exact description of how 3-D plasma shaping affects quantities 
such as the local magnetic shear and normal curvature which play crucial roles in the 
stability of localized plasma instabilities. In this work, this model is used to elucidate 
the mechanism by which small 3-D deformations can introduce substantial 3-D structure 
in the Pfirsch-Schluter current spectrum and local magnetic shear. 

A shaped tokamak equilibrium model, the widely used "Miller equilibrium" [26] . 
provides the starting point for these calculations. The physical position of the flux 
surface in cylindrical coordinates is parametrized in terms of the geometric poloidal 
angle 9 as 



where S is the plasma triangularity, k is the elongation, and A = Ro/r is the aspect ratio 
of the flux surface. The poloidal magnetic field on the surface is also specified using four 
additional parameters, 



where a = cos(xsin 9) + d r Ro cos 9 + [s K — s$ cos9 + (1 + s K )x cos#] sin 9 sm(9 + x sin 9), 
sin a; = 5 and three of the new parameters s K ,sg, and d r R$ are related to radial 
derivatives of the flux surface shaping parameters [26] . Typically q is chosen and d r ip is 
determined consistent with the definition q = (f /2ir) J dl p /(R 2 B p ) , where f(ip) = RB^, 
dl p is the differential poloidal arc length, and B^ is the toroidal magnetic field. The 
equilibrium parametrization is completed by choosing the following two dimensionless 
profile quantities, 



R = Ro + rcos[9 + (sin -1 5) sin 9], 
Z = nr sin 9 



(1) 
(2) 



d r i/j[sm 2 (9 + x sin + x cos 9) 2 + n 2 cos 2 9] 1/2 



(3) 
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where V' — J dQ j d^y/g/An 2 , yfg is the Jacobian and c =< B 2 Rq/\ViI)\ 2 > where 
< Q >= J J dQdC^Q^fgj J J dQdC,y/g. In total there are nine dimensionless parameters 
which uniquely specify a solution to the MHD equilibrium equations in the vicinity of 
the chosen flux surface. 

While the Miller equilibrium is formulated by denoting the inverse mapping X = 
X{9) and the poloidal field strength, the 3-D local equilibrium model is formulated in 
terms of straight field line coordinates and (, such that d(/dQ = q(ip). These two 
treatments can be unified by noting the relationship 
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Here, the Jacobian, yjg = R 2 / f, is written consistent with the use of symmetry in £. A 
transformation to the straight field line coordinates used by the 3-D local model can be 
calculated from 



dQ 



dX 



00 



(7) 



qR 2 B p 

The 3-D distortions of the flux surface shape can then be added perturbatively to 
the Miller equilibrium. In the following, the flux surface parametrization given by 

R = R(&) + Y,% cos(M,6 - Ni(), 

i 

Z = Z(Q) + sinlM.G - Ni(), (8) 

i 

and = — C will be used throughout the calculation, in which the sum is over 3-D 
perturbations with different helicity. Here, R(Q) and Z(Q) are given by Eqns ([1]) and 
and the amplitudes 7^ are small quantities. 

The magnitude of each 7, in the flux surface shape parametrization can be related 
to the magnitude of a radial magnetic perturbation. The resulting magnetic field is 
split into two components, B = ~Q axi -f B, where ~Q axi is the axisymmetric part. The 
perturbed part of the magnetic field is then projected into a radial coordinate, p 2 ~ ip axi . 
In a high aspect ratio circular cross section limit, the Fourier harmonic of each magnetic 
field component is given by 

% ^W-mA (9) 

-DO (M=Mi-l,N=Ni) n 

Note that at a rational surface q = M/N, the resonant component of the radial magnetic 
field is completely shielded. 

Using Eq. (j8]) as the flux surface parametrization, the geometric properties of 
the magnetic field lines (normal and geodesic curvature, normal torsion, etc.) are 
determined. The magnetic field and gradient in flux are given by 

1 f8X ldX\ . . 
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where the Jacobian is defined as 

ax /ax ax\ 
^ = a^'(ae x stJ ' (12) 

The local 3-D equilibrium model is based on an expansion in ip of the field line mapping, 

x(^, e, c) = x(^ , e, o + - ^0)^(^0, e, + ... (13) 

where ip labels the surface of interest and higher order terms in the expansion are 
neglected. By solving the ideal MHD equilibrium equations in the vicinity of ipo, dX/dip 
can be determined. To do this we note that the equations of MHD equilibrium dictate 
that no plasma currents flow normal to flux surfaces. By calculating n ■ J = 0, where 
n is the unit vector normal to the flux surface, a first order partial differential equation 
for the Jacobian is derived, 

d 9(( + 9ce/q _ d g & c + gee/q 



00 ^g d( ^g 



(14) 



where the metric elements are defined as gee = §H • §||) etc, and are completely defined 
by the parametrization given in Equation [2j Equation [14] then provides and equation for 
the Jacobian on the flux surface. The remaining components of dX/dip are determined 
by the conditions consistent with the MHD equililbrium conditions |23j . 



3. Ballooning stability 



The ideal ballooning stability boundary in s — a space for 3-D local equilibria can 
then be evaluated. These marginal stability boundaries can be constructed for 3-D 
configurations in a manner analogous to the traditional axisymmetric procedure [27] . 
The angular variables are transformed to another set which follows magnetic field lines, 
V = C (which marks the position along a magnetic field line) and «o — @ — C/q (the 
field line label). The lowest order incompressible ballooning eigenvalue equation at is 
written using R and V = Rq/B to normalize all physical quantities, 



d_ 
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Here, u 2 is a normalized eigenvalue, k 9 and n n are the geodesic and normal curvatures 
respectively (k = K n n + K g b x n), and A is the integrated local shear given by 
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\Vjj\ 2 V' 2 R 
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(16) 
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Figure 1. Spectrum of the radial magnetic perturbation for 3-D equilibria with several 
q values where N = 3 and ^ - 10 . 

Here, t/q plays the dual role of both the radial wave number of the mode as well as the 
starting point for field line integration. The quantity D is related to the local variation 
in the magnetic shear and will be treated in more detail subsequently One important 
feature of ballooning stability analysis in 3-D (relative to 2-D) is that the ballooning 
eigenvalue now depends on the field line label a®. In this work, separate marginal 
stability curves are calculated for different values of ocq. At each value of s and a, the 
ballooning eigenvalue equation is solved for values of rjk spanning one poloidal rotation of 
the magnetic surface to find the most unstable (or least stable) eigenvalue. In practice, 
for the 3-D equilibria considered here, the field line labeled by a® = has been found 
to be the most unstable field line on any given surface. 

The primary focus of this work is on the ideal ballooning stability properties of a 
family of equilibria calculated using this method. An axisymmetric equilibrium with 
A = 3.17, q = 3.03, k = 1.66, 5 = 0.416, s K = 0.7, s s = 1.37, and d r R Q = -0.354 is 
used as a base case. A set of equilibria are constructed with the same axisymmetric 
shaping, but with 3-D perturbations with 7^ = 0.001 , iVj = 3, and Mj = 4. .14 added. 
The different 3-D equilibria employ a series of q values approaching 3. Evaluating Eq. [9] 
for these parameters gives values of B p /Bq < 3 x 10 -3 for all harmonics of the equilibria 
examined here, with the M = 9, N = 3 component being of order 10 -5 (or even smaller 
as q approaches 3). The spectrum of B p / Bq is shown for several of the 3-D equilibria in 
Figure [TJ 

The marginal ballooning stability curves for an axisymmetric equilibrium as well 
as a family of 3-D equilibria with q ranging from 3.15 to 3.01 are shown in Figure 
Surprisingly, even very weak 3-D fields are capable of appreciably modifying the marginal 
stability boundary. The stability boundary is particularly sensitive to the q value, with 
the unstable region expanding as it approaches 3. This suggests that there is a resonant 
mechanism which affects the stability boundary, which is still weakly operative even 




Microinstability destabilization from shielded 3-D magnetic perturbations 



8 




Figure 2. Marginal ballooning stability curves for the axisymmetric equilibrium with 
q = 3.03 and equilibria with 3-D fields added for q values 3.15, 3.07, 3.05, 3.03, and 
3.01. All stability curves shown here are for the field line labeled by a = 0, which is 
generally the most unstable field line on each surface. 



when q = 3.15. It is also worth noting that the stability boundary expands dramatically 
at q = 3.01, which we will return to subsequently. 

The important question is then, how can a set of magnetic perturbations with a 
magnitude on the order of 10 -4 of the background field have such a dramatic impact 
on the ballooning stability boundary? We begin by inspecting all the terms appearing 
in Eq. ( fTBT) . A measure of how strongly various quantities are perturbed by the 3- 
D deformation is shown in Figure [31 The 3-D perturbation has no appreciable affect 
on B 2 ,\'Vip\ 2 ,Kg,K n or ^fg. However, the integrated local shear, A, is substantially 
modulated. By examining Eq. ( JT6l) it is clear that this modulation is due the local 
variation of the magnetic shear within the surface, encapsulated in the quantity D. 

The local variation of the magnetic shear is determined by a magnetic differential 
equation, 
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Figure 3. 

quantities. 

calculated: 



A measure of the degree to which the 3-D deformation modifies various 
For some quantity Q, evaluated on a real space grid, the following is 
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■ i.i 



1/2 



2D\ 
ij ) 



1/2 



A value of 0.01 corresponds to 



-a- 



1/2 

2q 2 



B 2 R 2 



A 



B 2 R 2 



B 2 Rl X 

|vv>| 2 v 7 



\v*Ij\ 2 v> |v# 



Co 



+2 



B 2 R 2 (^§r n Ro 



|V# 



b2r Kr 



c 



|V# 



(17) 



The right hand side of this expression can be split into three component parts relating to 
different physical effects. The first term on the right hand side relates to local shearing of 
the magnetic field due to changes in the average magnetic shear. The second term relates 
to shearing of the magnetic field due to Pfirsch-Schluter currents. The parallel current is 
written J\\/B =< J • B > / < B 2 > +p'\ where the quantity A is calculated consistent 
with the quasi- neutrality condition V • J = 0. The third term relates to shearing of 
the magnetic field due to the geometric properties of the surface (through the normal 
torsion r n = — n ■ (b ■ V)(6 x h)). The normal torsion is completely described by X(0, () 
[23] . Examining these three component parts shows that the 3-D perturbation does not 
appreciably affect the first term. There is a modest effect on the normal torsion but the 
dominant modulation of the local shear comes from the Pfirsch-Schluter currents. 

As a point of clarification, we note that the quantity D appears in the integrated 
local magnetic shear (Equation IT6]) . however it is dD/drj (the left hand side of Equation 
[T7|) which appears in the local magnetic shear. The resonant amplification of the right 
hand side of the Pfirsch-Schluter equation occurs once more when calculating D - this 
is due to the fact that the helical modulation of the local magnetic shear is nearly pitch 
resonant with magnetic field lines. Therefore, the integrated local magnetic shear (which 
explicitly appears in Equation [T5|) is modulated much more than the local magnetic 
shear. However, for the remainder of this work we will focus on the local magnetic 
shear, dD/dt], as it plays a more intuitive role in ballooning stability. 
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4. Local shear modulation 

The Pfirsch-Schliiter current spectrum is determined by a magnetic differential equation, 



Near rational magnetic surfaces, very small perturbations to the quantities on the 
right hand side of this equation can dominate the local magnetic shear. Using 
a Fourier decomposition, the Pfirsch-Schliiter harmonics are now given by A = 
Emat \mncos(MQ - N() with 



The differential operator in these magnetic differential equations becomes singular as q 
approaches a rational value, which can be seen as the denominator on the right hand 
side approaching zero. Even if the resonant component of the radial magnetic field is 
perfectly shielded, non-resonant components of the radial magnetic field can easily result 
in large responses to the right hand side of Eq. ( TT9l) . Figure H] shows the quantity D for 
several of the equilibria analyzed here. As q approaches 3, the peak amplitude of the 
variation of the local shear is greatly enhanced and the structure becomes dominantly 
helical. 

In order for ballooning modes to be driven unstable, the local magnetic shear must 
have a sufficiently low magnitude in regions of negative normal curvature. Figure [5] 
shows the distribution of the local magnetic shear for several 3-D equilibria with s = 1 
and a = 3, where the q = 3.01 equilibrium is near marginality but all other equilibria are 
well within the stable region. In the axisymmetric equilibrium, the local magnetic shear 
reaches relatively large negative values at the outboard midplane where the normal 
curvature is negative. The local magnetic shear passes through zero in a region of 
positive normal curvature. As q approaches 3, the helical modulation shifts the null 
in the local magnetic shear towards the outboard midplane. At q — 3.01, where the 
ballooning stability boundary is quite dramatically expanded, the null in the local shear 
overlaps with the region of negative normal curvature. This is strongly conducive to 
ballooning instability and explains the expansion of the unstable region in s — a space. 

If we limit our attention to the field line which passes through ( = and leave i]q 
fixed at 0, which is generally the most unstable point on a surface, more insight can be 
gained. First we examine the 3-D equilibrium with q = 3.05 at a = 4 and look at the 
effect of varying the surface average magnetic shear. Figure [6] shows the structure of 
the local magnetic shear as the flux surface averaged shear is raised from 2 to 5 (the 
stability boundary is at roughly 3.75), and the unstable eigenvector at the average shear 
of 4. As the surface averaged shear is raised, the strongly negative magnetic shear at 
the outboard midplane (G = 0) is weakened. The zeros in the local magnetic shear also 
move closer to the outboard midplane, where the unstable eigenvector peaks. 

This same process can be observed by examining equilibria with different q values 
at (s = 2,a = 4), where all equilibria except q = 3.01 are stable. Figure [7] shows the 






(19) 
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Figure 4. Contours of the quantity D on the magnetic surface for (a) the axisymmetric 
equilibrium with q = 3.03, and equilibria with 3-D fields added and q values of (b) 3.15 
and (c) 3.01. 

structure of the local magnetic shear for several equilibria, and the unstable eigenvector 
at q = 3.01. The local magnetic shear for an equilibrium with q = 3.005 is also shown 
(this equilibrium is so unstable that the stability boundary does not even fit in Figure 
EJ) As the safety factor approaches the resonant value, the local magnetic shear is 
helically modulated across the surface. For some magnetic field lines, this is catalytic 
for ballooning instability due to the lower magnitude of the local magnetic shear at the 
outboard midplane, and the shift in the zeros of the local magnetic shear towards the 




Figure 5. Contours of the local magnetic shear for (a) the axisymmetric base case, 
and 3-D equilibrium with (b) q = 3.07, (c) q = 3.03, and (d) q = 3.01. The thick black 
line marks the zero contour of the local magnetic shear, the thick blue line marks the 
zero contour of the normal curvature (which is negative near = 0), and the thick 
white line shows the path of a magnetic field line which passes through (0 = 0, C = 0). 

2 I ! ! ! ! ! 1 




Figure 6. 3-D equilibrium with q = 3.05, with a = 4. The local magnetic shear is 
shown at different values of the surface averaged shear. The ballooning eigenvector is 
also shown for the unstable s = 4 case. Here, the point of marginal stability is at a 
surface averaged shear of ~ 3.75. In all cases, ao = 0. 



outboard midplane. 
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Figure 7. The local magnetic shear for various equilibria with s = 2 , a = 4 , ao = 0. 
The unstable eigenvector is also shown for the q — 3.01 equilibrium. A 3-D equilibrium 
with q — 3.005 is also shown here to demonstrate how strong the local shear modulation 
becomes near resonance. 

However, for some of the magnetic field lines on the surface the result of the 3-D 
perturbations is to draw the zeros of the local magnetic shear away from the outboard 
midplane. Therefore the 3-D perturbations have a destabilizing effect on some areas of 
the surface but a stabilizing effect on others. Figure El shows the change in the structure 
of the local magnetic shear as the field line label is varied (770 is also varied such that 
the computational domain is centered at the outboard midplane). As the toroidal angle 
is varied away from ( = 0, the magnitude of the local magnetic shear is strengthened 
at the outboard midplane and the zeros shift towards the inboard side. A scan of the 
ballooning eigenvalue with respect to the toroidal angle is shown in Figure [9J showing 
that even deep within the unstable region, some portion of magnetic field lines are still 
ballooning stable. 

In principle, for 3-D configurations the stability of the entire flux surface should 
be studied self consistently. The infinite-n stability analysis is suitable in a tokamak 
where every magnetic field line on a surface has the same properties. However in 3- 
D geometry, as we have seen, different field lines can have radically different stability 
properties, so the analysis of a single field line may not be sufficient. In future work we 
will test the infinite-n MHD ballooning predictions with gyrokinetic KBM calculations 
which consider the full surface self-consistently. 

To conclude the stability analysis, we remark that the mechanism by which the 
local shear modulation affects ballooning stability is surprisingly simple. The helical 
perturbation to the local magnetic shear is nearly aligned with magnetic field lines, so 
an infinite-n ballooning mode located on a particular field line will mostly experience a 
uniform increase or decrease in the local magnetic shear (relative to 2-D) when the 3-D 
perturbation is added. This can have a stabilizing or destabilizing effect, for the exact 
same reasons that changing the surface averaged magnetic shear can be stabilizing or 
destabilizing depending on the location in s — a space. We have identified and elucidated 
this mechanism, however to perform quantitative predictions will require actual Kinetic 
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Figure 8. Change in local magnetic shear with field line label for the 3-D equilibrium 
with q = 3.01, a = 4, s = 2. The field lines labeled ao = and «o = 7r/15 are unstable 
while the other 3 are stable. The 3-D perturbation is destabilizing for some field lines 
but stabilizing for others. 
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Figure 9. The ballooning eigenvalue as a function of the field line label for the 3-D 
equilibrium with q = 3.01 at several values of the surface averaged shear. Here again, 
770 was varied to keep the computational domain centered at = 0. s = 1.5 is close to 
marginality and the s = 2,4 calculations go deeper into the unstable region. Even at 
s — 4, a portion of the magnetic field lines are still stable. 



Ballooning Mode calculations which consider an entire flux surface self-consistently, 
which is planned as future work. 

5. Summary and Discussion 

These calculations demonstrate two important points concerning the use of externally 
applied, 3-D RMP fields. The first is that very small 3-D magnetic perturbations, of 
experimentally relevant magnitudes, can drive substantial non-axisymmetric modulation 
of the Pfirsch-Schluter current spectrum and local magnetic shear even when shielding 
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physics prevents the presence of magnetic island formation. This effect is caused by a 
pitch resonance between magnetic field lines and the 3-D geometric deformation of the 
surface, in particular through a product of the geodesic curvature and magnetic field 
strength. This is related to the same physical effect which leads to the singular current 
problem of 3-D global MHD equilibrium calculations. 

A second point is that this 3-D modulation of the local shear has a significant 
effect on the ideal ballooning stability boundaries. These calculations suggest that 
the onset conditions for KBM instability can be significantly lowered near rational 
surfaces in the presence of 3-D perturbations. The initial variant of the predictive 
pedestal model EPED (EPED1) described the onset of KBM turbulence using the 
approxmation that the marginal stability boundary was of the shape OL crit ~ l/s 1 / 2 
[28] . suggesting that pedestal evolution is particularly sensitive to the lower left portion 
of the stability boundary, which is where it is most strongly perturbed by the 3-D 
local shear modulation. Gyrokinetic modeling of pedestal evolution has found that 
the inward advance of the pedestal occurs when unstable KBMs at the pedestal top are 
pushed toward marginal stability by the lowering of the surface averaged magnetic shear 
associated with growing boostrap currents as the pressure gradient rises [21]. Given that 
the dominant effect of the 3-D local shear modulation is to allow instability at much 
lower surface averaged shear, as seen in Figure [2J the 3-D pertubations could halt the 
inward advance of the pedestal by disrupting this stabilization process at the pedestal 
top once it reaches a low order rational surface. In this event, type-I ELMs can be 
avoided. 

There are several important consequences for future modeling work. Studies of 
plasma microinstabilities in the presence of 3-D RMP fields which utilize 2-D MHD 
equilibria and include the RMP field separately will not capture the Pfirsch-Schliiter 
physics. Proper analysis of this effect requires 3-D MHD equilibrium calculations with a 
finite pressure gradient at the surface(s) of interest. The ideal MHD model exhibits 
a singular parallel current at rational surfaces which in principle must be resolved 
using a more detailed physics model. We have sidestepped this by using radially local 
equilibrium calculations and merely examining 'near-resonant' effects. However, a better 
understanding of the nature of the plasma response to 3-D deformation at rational 
surfaces is desirable. 

In summary, we have presented a model for enhanced transport in the edge region 
of H-mode tokamaks in the presence of RMP fields. Small 3-D distortions of the flux 
surface shape can lead to helical Pfirsch-Schliiter currents which substantially modulate 
the local magnetic shear, a quantity which plays an important role in stability of localized 
modes and their associated anomalous transport. Even when the resonant components 
of the radial magnetic perturbation is shielded, this effect is still operative. The Pfirsch- 
Schliiter current and local shear are determined by magnetic differential equations which 
becomes singular near rational values of q. At a rational surface this singular response 
will be resolved by physical mechanisms beyond the scope of ideal MHD. However the 
near- resonant response is quite significant, and therefore 3-D magnetic perturbations 
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with a broad poloidal mode spectrum may be able to halt the inward evolution of 
the pedestal as it approaches a low order rational surface by adversely affecting KBM 
turbulence. 

Future work will use linear and nonlinear gyrokinetic calculation of KBMs to test 
the ideal MHD ballooning prediction. A better understanding of what physical processes 
resolve the singular Pfirsch-Schluter currents present in ideal MHD is desirable to 
assess the validity of the local 3-D MHD equilibrium calculations as a rational surface 
is approached. Sensitivity studies of these results with respect to external magnetic 
perturbation strength and spectrum content also need to be pursued. However, we have 
demonstrated that local magnetic shear modulation driven by helical Pfirsch-Schluter 
currents provides a possible explanation, which does not rely on stochasticity, for the 
observed enhanced plasma transport during some ELM suppression experiments. 



Appendix A. Analysis in the shifted circle limit 

In order to obtain gain more analytic insight into how the 3-D fields alter the ballooning 
stability properties, analysis of the ballooning equation for a particular 3-D equilibrium 
can be performed. We start from an axisymmetric equilibrium with large aspect ratio, 
low j3, circular concentric flux surfaces and locally steepened pressure gradient [29]. To 
this equilibrium, 3-D fields are added as described in Section II. Using Eqs. (1), (2) and 
(8) to define the equilibrium in the A 3> 1, 5 — k — s K — s$ — d r R = 0, ji/r <C 1 limit 
yields asymptotic estimates for various geometric quantities. In particular, the normal 
and geodesic curvatures are given by 

cos (~) *y- 

K n * S -——[l + O(-)] - £ M,-f- cosW - JViC), (Al) 
n i\ i n r 

^[1 + O(j)} + £ mn(M t - N t C). (A2) 



Kg 



Here, the first terms correspond to the usual 1/R contributions to the curvature vector, 
while the last terms correspond to the weak modulations due to the 3-D fields. While 
the 3-D fields provide very small corrections to the curvature vector, they can produce 
substantial corrections to the Pfirsch-Schluter coefficient. The solution to Eq. (18) for 
this equilibrium is given by 

A « 2^'^ [-goose - £ - cos(M t 9 - NtO], (A3) 

where the first term is the conventional high-aspect ratio prediction for the Pfirsch- 
Schluter coefficient and the last term is due to the 3-D fields. As noted, the 3-D fields 
can produce an appreciable effect when the 3-D field is near pitch resonance with the q 
profile. 

With these modifications, the 3-D corrections to the shifted circle equilibrium can 
be included assuming 7,/r 1 but (7j/r)(Mj — Niq)^ 1 ~ O(l). The ballooning equation 
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for the shifted-circle equilibrium with small 3-D distortions is given by 
8 BP 

— [(1 + A 2 )^|] + «[cos 6 + A sin 9]£ = -(1 + A 2 )^ 2 £, (A4) 

where the effect of the 3-D fields enters through the integrated local shear quantity A 
given by 

A= / d e{s-a(cose + ^-— -L_itcos[(M i -7V i g)e-JV i a )},(A5) 

where «o denotes a field line label. Generally, numerical solutions show that the large 
variation of the local shear produced by the 3-D fields produces localized eigenfunctions 
along the field line. Using this strong ballooning approximation, the ballooning equation 
can be converted to an equation of the form d 2 Y/dQ 2 + V(Q)Y = —u 2 Y, where the 
potential function V can be expanded about 6 = 0, V = V(0) + V"(0)Q 2 /2 + .... From 
this calculation, the corresponding critical a for ballooning instability can be calculated. 
In the large a, s limit, this is given approximately by 

a crit ~ — — ^ 777 V ± ( A6 ) 

l + EiM^- q ^cos(N iao ) 

where higher order corrections are suppressed. From this, we can see that the 3- 
D fields can produced order unity corrections to ballooning stability boundary when 
7i/r ~ (Mi — Niq). When this is the case, the a cr u can be substantially modified. 
Additionally, we note that the 3-D correction is sensitive to field-line label as noted by 
the ao dependence. Both of these features are also demonstrated by the more detailed 
numerical study described in the bulk of the paper. 
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